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Abstract. The coupling of the electromagnetic field to gravity is an age-old problem. 
Presently, there is a resurgence of interest in it, mainly for two reasons: (i) Experimental 
investigations are under way with ever increasing precision, be it in the laboratory or by 
observing outer space, (ii) One desires to test out alternatives to Einstein's gravitational 
theory, in particular those of a gauge-theoretical nature, like Einstein-Cartan theory 
or metric-affine gravity. — A clean discussion requires a reflection on the foundations 
of electrodynamics. If one bases electrodynamics on the conservation laws of electric 
charge and magnetic flux, one finds Maxwell's equations expressed in terms of the 
excitation H — {T>, Ti.) and the field strength F — {E, B) without any intervention of 
the metric or the linear connection of spacetime. In other words, there is still no coupling 
to gravity. Only the constitutive law H = functional(_F) mediates such a coupling. We 
discuss the different ways of how metric, nonmetricity, torsion, and curvature can come 
into play here. Along the way, we touch on non-local laws (Mashhoon), non-linear ones 
(Born-Infeld, Heisenberg-Euler, Plebaiiski), linear ones, including the Abelian axion 
(Ni), and find a method for deriving the metric from linear electrodynamics (Toupin, 
Schonberg). Finally, we discuss possible non-minimal coupling schemes. 

1 Introduction 

General relativity was proposed in 1915. One of its predictions was the bending 
of light rays of stars in the gravitational field of the Sun. This effect was verified 
observationally soon afterwards by Dyson et al. in 1920 and put, as a result, 
Einstein's theory in the forefront of gravitational research. 

Within the framework of general relativity, a light ray can be extracted from 
classical electrodynamics in its geometrical optics limit, i.e., for wavelengths 
much smaller than the local curvature radius of space. Accordingly, the bending 
of light can be understood as a result of a nontrivial refractive index of spacetime, 
see Skrotskii et al. ||8|,^ , due to the coupling of the electromagnetic field F to 
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the gravitational field g. Classically, we have in nature just these two fundamen- 
tal fields F and g, the weak and the strong fields being confined to microphysical 
dimensions of 10~^^m or 10~^^m, respectively. Therefore, the coupling of F and 
g is of foremost importance in classical physics. 

The conventional way that coupling is achieved is to display the Maxwell- 
Lorentz equations of vacuum electrodynamics in the (fiat) Minkowski world of 
special relativity theory in Cartesian coordinates. For this purpose, usually the 
formalism of tensor analysis (Ricci calculus) is used, see [f78| : 

F"' J = I" , Fv,k + Fjk,, + Fu^,j = . (1) 

Here Fij = -Fji = (i^oi, ■^02, ^bs, -F23, -F31, ■F'12) = {E,B) is the electromagnetic 
field strength, P the electric 4- vector current, and 

(2) 

with (7*^ as the contravariant components of the metric. The commas in 
denote partial differentiation with respect to the local spacetime coordinates x*. 

If we switch on gravity, the fiat Minkowski world becomes curved, the space- 
time geometry now being a Riemannian geometry with a variable metric gij {x) 

of Minkowskian signature (H ). The coupling of the Maxwell-Lorentz set 

(0) to gravity is now brought about by the comma goes to semicolon rule , — > ; 
(see [p2| ), where the semicolon represents the covariant derivative Vi = ;i with 
respect to the Riemannian connection ("Levi-Civita connection"): 

F'',j=I\ i^,;fc + + - . (3) 

This translation rule from special to general relativity is also alluded to as min- 
imal coupling with the additional understanding that the components of the 
metric in (|^) become spacetime dependent fields. 

The metric field gij(x), entering (|^) via (^ and via the covariant derivatives, 
i.e., via the semicolons, has to fulfill the Einstein field equation. 



Ricij 



1 h /Max mat \ ^ , 

-g^jRiCk'^ =n[^T ^, + T»jj , (4) 



Ric,j Rk^, , t\ ^ . ( - F.^F^'^ + \ Sf F^F'^' ] . (5) 



with 



MO V 4 

Here Riji'' is the curvature and T ij the material energy- momentum tensor. The 
coupled Einstein-Maxwell system describes correctly a wealth of experiments, in 
particular the gravitational bending of light, the gravitational redshift, the time 
delay of radar pulses in the gravitational field of the Sun, and the gravitational 
lensing and microlensing of starlight in the gravitational field of galaxies. 

But in all these experiments, we study the propagation of light along nuU- 
geodesics in a prescribed (and perhaps slowly varying) gravitational field which is 
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a solution of the Einstein vaeuum equation - and not of the electro-vacmxm equa- 
tion. We could call this the non self-consistent Einstein-Maxwell theory. In the 
solar system, e.g., the Schwarzschild metric is taken as solution of the Einstein 
vacuum equation and the motion of a "photon" is described by the null geodesic 
equation on this background. A true novel effect of the Einstein-Maxwell theory 
would be, e.g., the generation of electromagnetic waves by gravitational waves. 
Because of their smallness, no such effects were ever observed. Accordingly, the 
interaction of a classical electromagnetic field Fij in the form of a lightray with 
a prescribed gravitational field gij{x) is well described by means of Eqs.(§). 
Nevertheless, further consequences of these equations need to be compared with 
experiment as soon as more sensitive measuring methods are available. 

2 On the equivalence principle 

According to Einstein's equivalence principle, see [Q, gravity can be locally 
simulated in a gravity-free region of spacetime by going over from the Cartesian 
coordinates, anchored in an inertial frame of reference (including an intertial 
clock) and used in (|l|) , to arbitrary curvilinear coordinates yielding a non- inertial 
frame in general, as in (|^). In this context, the metric g^, occurring in and in 
the semicolons of (|[), is understood as a flat metric in curvilinear coordinates. 
Thus, the minimal coupling can be interpreted, in a first step, just as a coordi- 
nate transformation from Cartesian to curvilinear coordinates. And, moreover, 
it identifies the metric as the gravitational potential. 

On the other hand, let us assume that we are in a region with gravity and 
and ^ are valid together with the Einstein equation for the metric. Then, 
also according to Einstein's equivalence principle, we must be able to pick suit- 
able coordinates such that locally the equations look like in special relativity 
in Cartesian coordinates. In Riemannian geometry, the local coordinates are 
called Riemannian normal (hence geodesic) coordinates at one point P, if the 
Christoffcl symbols 

P^i'' := \ g'"^ {9u,3 + 9jLi - 9ij,i) (6) 

vanish at P and and the metric becomes Minkowskian: 

r,>|p^O, 5,,|p=diag(+l,-l,-l,-l). (7) 

Accordingly, the semicolon becomes a comma and the metric in (|^) , at one given 
point, looks flat. 

Still, the curvature is non-vanishing, of course: Rijk\p ^ 0. The equations 
look flat since they contain only first derivatives. If they contained second deriva- 
tives, then the semicolons goes to comma rule and its reverse would not work 
since on that level not only the Christoffels enter but potentially also the curva- 
ture which, in contrast to the Christoffels, is a tensor and cannot be nullified by 
means of a suitable choice of coordinates. For that reason, the minimal coupling 
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procedure, being in this context an expression of the equivalence principle, must 
be applied only to first order differential equations. The safest thing is then to 
apply it, as a rule, only on the level of a Lagrangian, since there ordinarily only 
first-order expressions are allowed for. Non-minimal couplings of the gravita- 
tional field to electromagnetism have also been investigated, see Prasanna fz^ , 
Buclidahl Goenner and Miiller-Hoissen | |5^ , |59| , |60[ , for example, or for 
light rays in non- minimally coupled theories, see Drummond and Hathrell p7[ |, 
but the price one has to pay is to introduce a new constant of nature; and there 
is no evidence for such a constant in nature - unless one takes the Planck length 
itself. We will come back to these questions in Sec. 7. 

Therefore we can conclude that the equivalence principle and minimal cou- 
pling work well for the Maxwell-Lorentz equations (0) and that they lead to 
experimentally established equations. 

Wave equation for the electromagnetic field strength We hasten to add 
that, within the framework of the minimally coupled Maxwell-Lorentz equations, we 
find 2nd derivatives if we derive the wave equation for the electromagnetic field strength 
F — and then also curvature terms are expected to emerge. This is exactly what 
happens, as already found by Gordon and Eddington |lq ]. 

In the framework of exterior calculus (see Frankel |2^), let us consider the elec- 
tromagnetic field strength 2- form F = ^ Fij dx^ A dx^ . In Maxwell-Lorentz vacuum 
electrodynamics, it satisfies dF = 0, eod *F = - J. We denote the codifFerential by 
S :— *d*. Then we find, with the wave operator (d'Alembertian) 

n~Sd + d5, (8) 

and by using the Maxwell-Lorentz equations, the wave equation 

nF= —d*J, (9) 

see [^. The left hand side of this equation, in terms of components, can be determined 
by substituting (|^: 

□ F = i (V" Vfc F,j + 2 Ric F,]fe - R'^.j F^i) dx' A dx' . (10) 

Accordingly, minimal coupling can lead to curvature terms of a prescribed form.^ 

3 A caveat 

Soon after general relativity had been proposed, it became clear, see Einstein 
[ p!9| , that one can introduce as auxiliary variables the densities 

^/^^9^\x)g^\x)Fu , X* := ^f^) P , (11) 



^ Of course, we could have non-minimal coupling as, e.g., in □F+7 (e^J e^gj i?"*^) AF 



— d *J , see Sec. 7. 
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with g(x) :— detgij{x), in terms of which the Maxwell-Lorentz equations (|^) can 
be rewritten in a metric-free way as 

:F'^ .j = I , F,J,k+F,k,^+ Ffc, J = . (12) 

Similarly, the charge conservation law 7*.^ = can be put in the form 

X\,=0. (13) 



The metric enters only via the densities defined in (|ll|). In fact, if we started from 
the set (12) in special relativity right away, then no comma goes to semicolon 
rule would have been necessary: These equations are generally covariant, they 
are valid in arbitrary curvilinear coordinates, be it in the framework of special 
or general relativity theory. 

In the calculus of exterior differential forms (Cartan calculus), see Frankel 
[ pT| , these equations can be formulated very succinctly. We introduce the electric 
current as odd 3-form, 

J p - j A = Jijk dx^ A dx^ A dx'' , (14) 
the electromagnetic excitation as odd 2-form 



H = V-n/\dt^ ^H.^dx' A dxi , (15) 



and the electromagnetic field strength as even 2-form 

F = B + E Adt = ^F.jdx' A dx^ . (16) 

Then (H) reads 

dH = J, dF = 0, (17) 

with 

dJ = 0. (18) 



The set (g/7|) represents the Maxwell equations. They are independent of metric 
and connection. The constitutive relation for the vacuum reads 

H=.[^*F, (19) 
V Mo 

where the star * represents the metric-dependent and odd Hodge duality oper- 
ator. Eq.(|9|) corresponds to (pl])i and I' can be related to the components of 
J, 

T'' = 4 \ , T = ^, e^-"^' JjM , (20) 

with e*-''^' = ±1,0, the totally antisymmetric Levi-Civita tensor density. 

Now the equivalence principle looks empty: Since the Maxwell equations ( p^ 
are formulated in a coordinate and frame independent way, they are valid in this 
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form in arbitrary coordinate systems and frames, be it in a flat or in a curved 
spacetime. Only the constitutive relation (|l9|) "feels" , up to a conformal factor, 
the presence of a flat or a non-flat metric, i.e., the constitutive relation couples to 
the conformally invariant part of the metric. The coupling of electromagnetism 
to gravity becomes almost trivial. Is all this just a mathematical trick, which 
distracts from the physical content of Maxwell's theory, or is it more? 

One further observation hints also at the need for clarification. The Einstein- 
Cartan theory of gravity is a viable gravitational theory, see |3^ , |5l|,^. It is the 
simplest model of the metric-affine gauge theory of gravity, see |3l| , |25|] . In the 
Einstein-Cartan theory, spacetime is described by means of a Riemann-Cartan 
geometry with torsion and curvature]^ If we couple (|^) to gravity, do we have to 
use the semicolons as covariant derivatives with respect to the Riemann-Cartan 
connection or still with respect to the Christoffels, see [0? In the context of 
(^^, this question cannot even be posed, since the exterior derivative d is all 
what is needed. Are then the equations (|l^) misleading as a starting point for 
coupling to gauge gravity? What could be the appropriate starting point? 

Provided one formulates Maxwell's theory and its coupling to gravity in 
terms of a Lagrangian with the electromagnetic potential A as variable, gauge 
covariance of the formalism results in ( [l7| ) cum ([l9|), as was pointed out by Benn, 
Dereli, and Tucker However, we would like to have some more immediate 
insight into the structure of electromagnetism as induced by experiment even 
without having a variational formulation at our disposal. 



4 Electric charge and magnetic flux conservation 



The metric is a quantity which allows to define lengths and angles in spacetime. 
There are, however, laws in physics which don't require the knowledge of a 
metric. Take the conservation law of electric charge as an example. Mark a 3- 
dimensional simply connected submanifold ^2^. We know from experiment that a 
possible electric charge inside is composed of charge "quanta", i.e., there is an 
integer number of elementary charges in fl^^. Recent advances in technology made 
it possible, see [|l3[|4|, to trap and to count single electrons and protons. Thus, 
as soon as we have such quanta available, we can rely on counting procedures, 
see Post l?^ , the use of a meter stick or a chronometer is superfiuous under such 
circumstances. 

Electric charge conservations is experimentally well-established and is one of 
the pillars electromagnetism rests on. We formulate it, following Kottler-Cartan- 



van Dantzig, see |72,|87[ and also |32|, most appropriately as an integral law. 
According to we assume the existence of the odd electric current 3-form J. 
We take charge conservation as axiom 1, that is, J integrated over a closed 3- 
dimensional hypersurface fl^ has to vanish, if this hypersurface is the boundary 

^ A proper discussion of the equivalence principle in the context of Einstein-Cartan the- 
ory requires the introduction of local coframes, see [ p5| , p7p^ ] . Being concerned here 
only with electromagnetism, it is sufficient to use natural, i.e., holonomic coframes. 
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of a connected 4- volume f24: 



J ^ J dJ ^0, or dJ ^0. (21) 

Here we applied the Stokes theorem. 

If (^ ) is assumed to be valid /^^ J = for all three-cycles 9C3 = 0, C3 ^ 



dfli, then, according to a theorem of de Rham, J is exact, see |73 . Thus the 
inhoniogeneous Maxwell equation is a consequence, 

J = dH , (22) 

with the odd 2-form H of the electromagnetic excitation, see (|l5| ) . The excitation 
is only determined up to an exact form. Nevertheless, the electric excitation T) 
can be measured by means of Maxwellian double plates as charge per unit area, 
the magnetic excitation Ti. by means of a small test coil, which compensates 
the Ti-field to be measured, as current per unit length. This is possible since in 
theses null experiments vanishing field strength F implies vanishing excitation 
H, see |3^. In other words, the extensive quantities V and Ti. - and thus the 
4-dimensional excitation H - have an operationally significance of their own, 
since they are related to charge at rest or in motion, respectively. Accordingly, 



the somewhat formalistic introduction of the densities in (11) has now been 
legitimized as a transition to operationally meaningful additive quantities. Note 
that up to now only the differential structure of the spacetime was needed, a 
metric has not been involved. 

Let us choose a field of 4- frames Ca and consider the motion of a point particle 
with respect to the reference frame thus defined. As axiom 2 one can take an 
operational definition of the electromagnetic field strength F via the Lorentz 
force density 

U = iec.\F)AJ. (23) 

The interior product (contraction) is denoted by J . The force density fa is a 
notion from classical mechanics. It is an odd covector- valued 4-form. Accordingly, 



Eq.(|23|) can be read as a definition of the even 2-form F, see (|16D. Again, we 
don't need a metric. And we know the recipe of how to proceed in the same 
manner. 

That F can be interpreted as magnetic flux is obvious if we choose f22 as 
a 'spacelike' surface (strictly, at this point we don't know what spacelike means; 
we will come back to this later). In superconductors under suitable circumstances 
we can count (in an Abrikosov flux line lattice) quantized magnetic flux lines. 
This suggests that magnetic flux is a conserved quantity (axiom 3): 

F = J dF ^0 or dF ^0. (24) 



In this way, by means of the axioms (|21|), (|23|), and (^4|), we recovered the 
fundamental structure of Maxwell's theory: dH — J, dF = . This is what had 
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been called metric-free electrodynamics^ What is missing so far is the relation 
between the excitation H and the field strength F, and it is exactly there where 
the metric, i.e., the gravitational potential comes in. 



5 No interaction of charge and flux 
"substrata" with gravity 

We now understand that the inhomogeneous Maxwell equation dH ~ J, as 
an expression of electric charge conservation, cannot be influenced by gravity, 
i.e., by the metric tensor g, or, in the case of metric-affine gravity or its specific 
subcases, such as Einstein-Cartan theory, by the connection F of spacetime. The 
electric charge substratum of spacetime has rules of its own. Spacetime can be 
deformed by the presence of metric and connection, but the charge substratum 
and the net electric charge to be attributed to a prescribed 3-dimensional (3D) 
volume won't change. Thus the additivity 3D volume-wise of the charge lays 
at the foundation of the Maxwellian framework. And it translates into the 2D 
additivity of the integrated excitation J^^ H - this being the reason why one 
uses this integral for the operational interpretation of H. 

Similar arguments can be advanced for the homogeneous Maxwell equation 
dF — 0. However, first of all it should be stressed that the axiomatics we are 
using strongly suggests the non-existence of magnetic charges. If there were 
magnetic charges, then we would have no reason to believe in electric charge 
conservation either; compare for this argument axiom 1, Eq.(|l|), with axiom 3, 



Eq.(24). Conventionally, the inhomogeneous equation dH = J is seen in analogy 
to the homogeneous one dF = 0. But not so in the framework of our axiomatics 
which has a firm empirical basis. We put d J = in analogy to dF = 0. The whole 
historical development of electromagnetism, starting with 0rsted and Ampere, 
points to the elimination of the phenomenologically introduced magnetic charges. 
Most recent experiments, see ]i|,p9t, exclude magnetic charges with very good 
precision. Thus theoretical as well as experimental evidence speak against the 
existence of magnetic charges. 

Having said this, we hasten to add that, nevertheless, there is some kind 
of magnetic substratum in spacetime, namely the magnetic flux J^^^ F. It is a 
substratum of its own right. The fluxoids, the quantized magnetic flux lines in 
superconductors, see [Q, do convey a clear message. Besides electric charge^ 
magnetic flux^ (and not magnetic charge) has an independent standing in elec- 

^ Stachel |Q calls it generalized electrodynamics. We don't follow this suggestion, 
since Maxwell's equations were originally given in terms of {T>, TC) and {E, fiH) in 
a form 'isomorphic' to the (1 -I- 3) -decomposition of dH = J and dF = 0, see |5C[|. 
Therefore, the "generalized" Maxwell equations, dH — J and dF = 0, correspond 
in actual fact, just to Maxwell's equations (modulo the substitution B fiH). And 
this is how we will name them. 

Sl-unit Coulomb, elementary charge e = 1.60217733 x 10"^'' C. 
^ Sl-unit Weber, elementary fluxoid h/{2e) = 2.06783461 x 10"^^ Wb. 
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tromagnetism, too. Thus rightfully, it is governed by an own axiom, namely 
axiom 3. 

Axiom 3 is again a conservation theorem. In contrast to axiom 1, which has 
a fermionic smell, axiom 3 is more of a bosonic nature. Moreover, magnetic flux 
adds up 2D area-wise. For this reason, magnetic flux is represented by a 2-form 
and not, like the charge, by a 3-form. Accordingly, there are essential differ- 
ences between these two conservation laws which express the peculiarities of the 
electromagnetic phenomena. Electric and magnetic effects enter the Maxwellian 
framework in an asymmetric way, in spite of all that talk about a duality between 
electricity and magnetism. But there is also a similarity in that both axioms are 
formulated as integral conservation laws. The possibility to count the fluxoids 
assures us that axiom 3 has to be again a law free of metric and connection. 

Incidentally, there is a nice visualization of the fundamental quantities en- 
tering electrodynamics. If one describes the quantum Hall effect for low lying 
Landau levels, then the concept of a composite fermion is very helpful: it consists 
of one electron and an even number of fluxoids is attached to it, see Jain pj],[42| 



and 1 43 . Isn't that a very clear indication of what the fundamental quantities 
are in electrodynamics? Namely, electric charge (see axiom 1) and magnetic flux 
(see axiom 3), see also Nambu in this general context. 

Our conclusion is then that, as long as we opt for electric charge and mag- 
netic flux conservation, the Maxwell equations in gravity- free regions, i.e., in 
the Minkowski spacetime of special relativity, read dH ~ J and dF ~ 0; they 
remain the same irrespective of the switching on of gravity, be it in Einstein's 
theory, in metric-affine gravity (see |]76|), or in any other geometrical theory of 
gravity. 



6 Constitutive law of electrodynamics and its relation to 
gravity 

After having discussed extensively that gravity does not influence the Max- 
well equations, we eventually turn to the constitutive law via which gravity 
does influence electrodynamics. It is true, the charge substratum and the flux 
substratum themselves do not couple to gravity, as we have shown in the last 
section. However, the interrelationship between both substrata is affected by 
gravity. Metaphorically speaking, the "flow" of each of the substrata is ruled by 
a particular gravity-independent conservation theorem, but the flows of electric 
charge and magnetic flux are coupled via a gravity-dependent constitutive law 
since, in the end, magnetism has to be expressed in terms of electricity. 
Let us choose arbitrary local spacetime coordinates x^. Then we have, 

= i iJy dx' hdx\ F = i dx' A dx^ . (25) 

We will turn flrst to the electrodynamics of material media in order to develop 
some intuition on the concept involved, but eventually, it will be the vacuum^ be 
it in inertial or non-incrtial frames, which will occupy the center of our interest. 
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6.1 Non-local 

Moving macroscopic matter defines a (1 + 3)-splitting of spacetime specified by 
a well defined average 4-vector velocity field u which describes the congruence 
of worldlines of the flow of the medium. Such a vector field can be defined 
operationally from the motion of matter as follows. 

Let a 3-dimensional arithmetic space be equipped with the coordinates 
a = 1, 2, 3. We will use these coordinates (known as Lagrange coordinates in 
continuum mechanics) as labels which enumerate elements of a material medium. 
A smooth mapping xjq) : B? into the spacetime defines a 3-dimensional 

space domain (hypersurface) V which represents the initial distribution of mat- 
ter. In local spacetime coordinates, this mapping (or labeling) is given by the four 
functions xLj(^"). It should be preserved at any time, i.e. along any worldline 
of a particular element its labels are constant. 

Given the initial configuration V of matter, we parameterize dynamics of 
the medium by the "time" coordinate r which is defined as the proper time 
measured along an element's worldline from the original hypersurface V . The 
resulting local coordinates (t, ^°) are usually called the normalized comoving 
coordinates. Thus finally, the motion of matter is described by the functions 
x^{t,^°-). Subsequently, we define the 4- velocity vector field by 



Evidently, a family of observers comoving with the matter is characterized by 
the same timelike congruence x*(r, ^°). They are making physical (in particular, 
electrodynamical) measurements in their local reference frames which drift with 
the material motion. 

One says that a medium, moving in general, has dispersion properties when 
the electromagnetic fields produce non-instantaneous polarization and magne- 
tization effects. The most general linear constitutive law is then given, in the 
comoving system, by means of the integral 



The coefficients of the kernel Kij^^ir^r') are called the response functions. We 
expect the metric to be involved in their set-up. Their form is defined by the 
internal properties of matter and by the motion of a medium. 

Mashhoon |Q has proposed a physically very interesting example of such a 
non-local electrodynamics in which non-locality comes as a direct consequence of 
the non-inertial dynamics of observers. In this case, instead of (^), one should 
use the field expansions 




—const 



(26) 




(27) 



F = 



(28) 
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with respect to the coframe of a non-inertial observer — Ci" dx^. The consti- 
tutive law is then replaced by 

Hap{r,0 dr'K^f,''\r,T')F^s{r',0 , (29) 

and the response kernel in (p9| ) is now defined by the acceleration and rotation 
of the observer's reference system. It is a constitutive law for the vacuum as 
viewed from a non-inertial frame of reference. 

Mashhoon imposes an addititional assumption that the kernel is of convolu- 
tion type, i.e., Kap^^ {t,t') = Ka/s'^^ir — r'). Then the kernel can be uniquely 
determined by means of the Volterra technique, and often it is possible to use 
the Laplace transformation in order to simplify the computations. Unfortunately, 
although Mashhoon's kernel is always calculable in principle, in actual practice 
one normally cannot obtain K explicitly in terms of the observer's acceleration 
and rotation. 

Preserving the main ideas of Mashhoon's approach, one can abandon the 
convolution condition. Then the general form of the kernel can be worked out 
explicitly (u is the observer's 4- velocity): 

i^a/'(r,r') = 1 e^0'^'{s-J^ <5(r - r') - u\n''\r')) . (30) 

The influence of non-inertiality is manifest in the presence of the connection 
1-form. The kernel (^0|) coincides with the original Mashhoon kernel in the case 
of constant acceleration and rotation, but in general the two kernels are differ- 
ent Perhaps, only the direct observations would establish the true form of 
the non-local constitutive law. However, such a non-local effect has not been 
confirmed experimentally as yet. 



6.2 Non-linear 

But the constitutive law can also be non-linear (or non-local and non-linear at 
the same time) . In the local and non- linear Born-Infeld electrodynamics |^ , with 
the dimensionfuU parameter /c as maximal attainable electric field strength, we 
have 

The metric as symmetric second rank tensor enters here in a very natural way. It 
adds up with the antisymmetric electromagnetic field to an asymmetric tensor 
- much in the way Einstein had hoped to find for his unified field theories of 
gravity and electromagnetism. By differentiation, we find 

*F + ^"{F hF)F 
H=J- , ' ' (32) 

-j,*{FA*F)-^[*(FAFW 
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now the metric being absorbed in the (odd) Hodge star operator, see For 
/c — > cxD, we recover the conventional local and linear Maxwell- Lorentz theory 
for vacuum with H — ^J^*F. The Born-Infeld electrodynamics is presently 

used as a toy model in string theories, see |Q. The problem with Born-Infeld 
electrodynamics is that, in contrast to Maxwell's theory, it defies quantization. 
It is an interesting model, but nothing like an established theory. 

A similar example is the non-linear Heisenherg-Euler electrodynamics ||3^ . 
Quantum electrodynamical vacuum fluctuations yield corrections to Maxwell's 
theory that can be accounted for by an effective constitutive law constructed 
by Heisenberg and Euler. To the first order in the fine structure constant af = 



' ' '" 8 Q!f 



, it is given by (see also |39[|35|) 

F h*F) 



45 Bl 



14 Off 



where Bk — ^^^^ ~ 4.4 x 10^ T, with the mass of the electron m. The met- 
ric is again hidden in the Hodge star and the Maxwell-Lorentz limit results 
analogously for to ^ oo. The Casimir force between two uncharged electrically 
conducting plates, also an effect of vacuum fluctuations, has been experimen- 
tally verified as have been non-linear effects in the "superposition" of strong 
laser beams. Accordingly, the non-linear constitutive law (^3|) is a valid post- 
classical approximation of vacuum electrodynamics and as such experimentally 
confirmed. 

Note that these variants of classical electrodynamics respect charge and flux 
conservation. This underlines the fact that our axiomatics clearly points to that 
structure of electrodynamics, namely the constitutive law, which can be changed 
without giving up the essentials of electrodynamics. 

Both, Eqs.(^2|) and (|3^ ) are special cases of Plebahski's more general non- 
linear electrodynamics [|7l| . Let the quadratic invariants of the electromagnetic 
field strength be denoted by 

h ■.^^*{F A*F) = ^{E^ - B^) and h := ^* {F /\ F) ^ E ■ B , (34) 

where /i is an even and I2 is an odd scalar (the Hodge operator is odd). Then 
Plebaiiski postulated a non-linear electrodynamics with the constitutive law^ 

H = U{h,l2)*F + V{h,h)F, (35) 

where U and V are functions of the two invariants. Note that in the Born-Infeld 
case U and V depend on both invariants whereas in the Heisenberg-Euler case 
we have C/he(^i) and Vhe(^2)- Nevertheless, in both cases U is required as well 
as V. And in both cases, see ( |3^ ) and (|33|), U is an even function and V and 
odd one such as to preserve parity invariance. 

® Strictly, Plebariski assumed a Lagrangian which yields (^^ together with the struc- 
tural relations F = u{h.l2)*H + v{Ii, I2) H . The latter law, apart from singular 
cases, is equivalent to (651). 
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If one chose V to be an even function, e.g., then parity violating terms would 
emerge. Such terms were most recently discussed by Majumdar, Mukhopad- 
hyaya, and SenGupta [ ^SUs^ ; for the experimental situation (there seem no 
signatures for parity violations) compare Lue et al. 

Singularity-free electro-gravitodynamics Recently, Ayon-Beato & Garcia 
[H , for earlier work see Shikin ||7^ , have proposed a constitutive law 

H = U{h)*F, (36) 



which, as subcases, does neither encompass ( p2D nor (|33| ) and thus makes it 
appear as rather academic. The explicit form of U is defined by the requirement 
of obtaining completely singularity-free solutions of the coupled system of the 
gravitational field (Einstein) and the electromagnetic field (non- linear Maxwell). 
Examples of suitable functions U{Ii) are given in (t^,^. 

In terms of the local time and space coordinates {t, r, 6, (p), the general spher- 
ically symmetric ansatz for the coframe can be written as 

i)^=f(r)dt, iji^-i-dr, ^^^rde, -d^ = r sm0 d(j> , (37) 
fir) 

whereas, for the electromagnetic field, we have 

F = ip{r) A^K (38) 

The exact solution of the coupled system of gravitational and electromagnetic 
field equations, i.e., of Einstein's equation (||) and Maxwell's equations dF = 
0, dH = 0, reads 

,n g f2 2m Q(r) 

(y(ii)r^ r 

where 9,™ are integration constants and ('Tolman's integral') 

oc li 

Q{r) = Kr [ dr' IC{r')r'^, K = 2IiU{h) - ( dI[U{I[). (40) 



In the last function one should substitute the explicit form of the quadratic 
invariant Ii computed on the spherically symmetric configuration ( p7| ) with (|38|). 

It is shown in |7^,||,||,^ that the constitutive function C/(/i) can be chosen 
in such a way that the functions in (^9|) describe a completely regular, i.e., 
singularity-free configuration. 

6.3 Linear: Abelian axion, inter alia 

A very important case is that of a linear constitutive law between the compo- 
nents of the two-forms H and F. It postulates the existence of the 6 x 6 = 36 
constitutive functions Hij^^{t, x) ~ — nji^^ = — Ki/'^ such that 

^^K,,''' Fki. (41) 
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This kind of an ansatz we know from the physics of anisotropic crystals. The 
factor 1/2 is chosen in order to have a smooth transition to the conventional 
D = SqeE etc. relations, cf. ||7|l p. 127. 

The choice of the local coordinates is clearly unimportant. In a different coor- 
dinate system the linear constitutive law preserves its form due to the tensorial 
transformation properties of Kij^K Alternatively, instead of the local coordinates, 
one may choose an anholonomic frame and may then decompose the two-forms 
H and F with respect to it. 

Since H is an odd and F an even form, the constitutive functions Kij'^^t, x) 
are odd. Taking the Levi-Civita symbol, we can split off the odd piece according 
to 

^..fci_. 1,.. mnkl ^ijkl _ }_ ijmn kl / Ar^\ 

Because of the corresponding properties of the Levi-Civita symbol, the x^''^^ £^re 
even scalar densities of weight -1-1. For the Levi-Civita symbols with upper and 
lower indices, we have e'^*^' emnpq = ^m^pq- 

With the linear constitutive law (as with more general laws), we can set 
up a Lagrangian 4-form; here we call it Viin. Because oi H = — dVun/dF, the 
Lagrangian must be quadratic in F. Thus we find 



'-H.F^-l 
2 8 



Viin = --H /\F = -- HijFpq dx' A dx' A dxP A dx« 



= - ^ (e^™ey>n„x""'') FuFp, dx° A dx' A dx^ A dx^ . (43) 

The components of the field strength F enter in a symmetric way. Therefore, 
without loss of generality, we can impose the symmetry condition x'"''^' — X^^^'' 
on the constitutive functions reducing them to 21 independent functions at this 
stage. 

The Kij^^ carry the dimension [k] — [x] — /fi. Therefore, still before intro- 
ducing the metric, we can split off the totally antisymmetric part of x^''^^ 
define the dimensionless constitutive functions according to 

X'^^^ = f X '^^^ + a e'^^^ , with X ^'^'^'^ = . (44) 

Here [/] = [a] — fi/e^, and / = fit,x) and a — a(t,x) represent one scalar and 
one pseudo-scalar constitutive function, respectively. Thus the linearity ansatz 
eventually reads 

u ^ mnkl TP f mnkl , jji /At:\ 

J^ij — ^^ijmnX i'kl — ^^ijmn X f M + OL t ij , ^4t)j 

with 

^ninkl ^nmkl ^mnlk ^klmn and ^[mnkl] q (46) 

o , , 

i.e., besides a, we have 20 independent constitutive functions. Thus x"™ has 
the same algebraic symmetries and the same number of independent components 
as a curvature tensor in a Riemannian spacetime. 
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Pseudo-scalars are also called axial scalars. So far, our axial scalar a{x) is 
some kind of permittivity/permeability field. If one adds a kinetic term of the a- 
field to the electromagnetic Lagrangian (|4|), then a{x) becomes propagating and 
one can call it legitimately an Abelian[^ axion. Ni |^ was the first to introduce 
such an axion field a in the context of the coupling of electromagnetism to 
gravity, see also deSabbata & Sivaram and the references given there. 

The Abelian axion has the following properties: 

• Pseudoscalar field, i.e., spin ~ 0, parity = — 1. 

• Couples to Maxwell's field in the Lagrangian according to aF/\F = 2a EA 
BAdt, see ( ^5|) and (|4^). Here E is the 3-dimensional electric field 1-form and 
B the corresponding magnetic field 2-form. This term in the total Lagrangian 
can be written as 

aFAF = -daAAAF, (47) 

dropping, as usual, the total derivative. This contributes to the excitation 
H = -dL/dF a term 

^ daAA. (48) 

• Since it arises on the same level as the metric, see Eq.(|^ below, it is a field 
of a similar universality as the gravitational field. 

As yet, the Abelian axion has not been found experimentally, see the discus- 
sion of Cooper & Stedman JT^ on corresponding ring laser experiments. 



6.4 Isotropic 

The linearity ansatz ( p5[ ) can be further constrained in order to arrive eventually 
at an isotropic constitutive tensor. We will proceed here somewhat unconven- 
tional in that we don't assume a metric of spacetime beforehand but rather 
derive it in the following way: 



Duality operator, electric and magnetic reciprocity The constitutive 

tensor X '''™" of ( ^ ) defines a new duality operator which acts on 2-forms on 
X. In components, an arbitrary 2-form = \Oij dx^ A dx^ is mapped into the 
2-form by 

see No metric is involved in this process. Now the linear material law 

( ^5| ) can be written as 

H = {f* + a) F. (50) 

We postulate that the duality operator, applied twice, should, up to a sign, 
lead back to the identity. Such a closure relation or the "electric and magnetic 

In contrast to the axions related to non- Abelian gauge theories, see | pl| , p2||5^ and 
the reviews in ^ and Q. 
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reciprocity" [ psj reduces the number of independent components of X to 9 (with- 
out using a metric). One can demonstrate that this is a sufficient condition 
for the nonexistence of birefringence in vacuum, see |6^ , |63|j4^ , |6^ , ^ . Then the 
fourth-order general Fresnel equation degenerates to the second-order hght cone 
equation. Therefore, we impose 

## = -1. (51) 

The minus sign yields Minkowskian signature^, whereas the condition "^"^ = +1 
would lead to Euclidean or to the mixed signature (-I-, — , — ). 

Seemingly Toupin and Schonberg j?^ were the first to deduce the confor- 
mally invariant part of a spacetime metric from duality operators and relations 
like (IH) and (|l]). This was later rediscovered by Jadczyk |^^, whereas Wang 
[ pO| gave a revised presentation of Toupin's results. A forerunner was Peres pOt , 
see in this context also the more recent papers by Piron and Moore ||7^ . Brans 
and subsequently numerous authors discussed such structures in the frame- 
work of general relativity theory, see, e.g., ||,^,^,^ and the references given 
there. 

It is convenient to adopt a more compact bivector notation by defining the 
indices I, J,... = (01,02,03,23,31,12). Then X^^*^' becomes the 6 x 6 matrix 
X and ( |5l| ) goes over into 

x"ejKX''''eLM^-Sii. (52) 

In terms of 3 x 3-constituents an arbitrary symmetric X =X constitutive 
matrix reads 



where A — A ,B = B , and the superscript denotes transposition. The 
° r 7 

algebraic condition e/j X = is provided by trC = 0. 

The general non-trivial solution of the closure relation (^) can be written in 
the form 

J" ^(^-;-^" -;*■). (M, 

Here i? is a nondegenerate arbitrary symmetric 3x3 matrix (6 independent com- 
ponents Bab), K an arbitrary antisymmetric matrix (3 independent components 
Kab ='■ ^abck'^), N the symmetric matrix with components N°-^ := k'^k'', and 
q := -1/det B, p := [tr{NB)/det B] - 1. Thus, Eq.(||) subsumes 9 independent 
components. 



One could define a different duality operator by "^Oij = jeijki x*'™" Omn such that 
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Triplet of self-dual 2-forms The duality operator induces a decomposition 
of the 6-dimensional space of 2-forms into two 3-dimensional invariant subspaces 
corresponding to the eigenvalues ±i. Writing the 2-form basis 0^ — dx* A dx^ in 
terms of the two 3-dimensional column vectors 

a,6,... = 1,2,3, (55) 

one can construct the corresponding self-dual basis ^{'d^ ~ * "^O^). In the 

3-vector representation, 

§'=(1), (56) 

one of the 3-dimensional invariant subspaces can be spanned either by the up- 
per or by the lower components which are related to each other by a linear 

(^) _ (s) 

transformation. For example, /? can be expressed in terms of 7 according to 

(s) _ _ (s) (s) 

(3 — (i + B K)B 7 . Therefore 7 or, equivalently, the triplet of 2-forms 

= ^ {dx° A dx" - (det B)'^ h dx^ A dx" 

+ i{B^^Y^ ebcddx" Adx"^). (57) 

subsume the properties of this invariant subspace. Each of the 2-forms carry 3 
independent components, i.e., they add up to 9 components. 

o 

The information of the constitutive matrix X is now encoded into the 
triplet of 2-forms 3^°"^ . One can verify that the latter satisfies the completeness 
relation 

S-^"^) A S*'^) = I {B-^y^ {B)cd S^"^ A S^'^^ . (58) 



Extracting the metric Within the context of SU{2) Yang-Mills theory, Ur- 
bantke was able to derive a 4-dimensional spacetime metric Qij from a triplet 
of 2-forms satisfying a completeness condition of the type (|5^). Explicitly, the 
Urbantke formulas read 

v/d^ g., =--\ Vd^ ea,c 6'=""" S'^S'^S';:^ , (59) 

x/d^ e''™" Sifsi;il . (60) 

The slj^ are the components of the 2-form triplet S^"-^ = s\f'dx^ A dx"^ /2. If we 
substitute the forms (^7|) into and (|6^) , we can display the metric explicitly 
in terms of the constitutive coefhcients: 
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Here ka := Babk^ = Bab^^^'^'Kcd/'^- One can verify that the metric in (61) has 
Minkowskian signature. Since the triplet 5'*^°-' is defined up to an arbitrary scalar 
factor, we obtain a conformal class of metrics. 

Given a metric, we can now define eventually the notion of local isotropy. 
Let T*i "*p be the contravariant coordinate components of a tensor field and 
rpai...ap eij"i ■ ■ ■ Bi "pT*i- -'p its frame components with respect to an or- 
thonormal frame Gq, = e*Q. di. A tensor is said to be locally isotropic at a given 
point, if its frame components are invariant under a Lorentz rotation of the 
orthonormal frame. Similar considerations extend to tensor densities. 

There are only two geometrical objects which are numerically invariant under 
local Lorentz transformations: the Minkowski metric Oap — diag(+l, —1, —1, —1) 
and the Levi-Civita tensor density eap-ys- Thus 

T'^'' = <Pix) {g^'g^' - g^'^g^') + ip{x) e«^'=' (62) 

is the most general locally isotropic contravariant fourth rank tensor density of 
weight +1 with the symmetries T*-"^' = —q'j'-ki _ _q-ijik _ rj-kUj _ jjgj.g ^ ^ 

are scalar and pseudo-scalar fields, respectively. 

One can prove that the constitutive tensor in fiS) with the closure property 
is locally isotropic with respect to the metric (^ip , see also . Accordingly, 
for the constitutive tensor, we finally have 

X*^'^' =2V^5'=V''- (63) 

Thus, the isotropic law reads 

H^, = ^ g'^™^'" Fki + a F,, (64) 

or, if written with the help of the Hodge star operator belonging to the metric 

H^{f* + a)F. (65) 

6.5 Centrosymmetric 

If we want the constitutive tensor to be refiection symmetric at each point of 
spacetime, i.e., if we require centrosymmetry, then we have to kill the Abelian 
axion and arrive, provided F is chosen in accordance with the Sl-conventions, 
at the usual law for Maxwell-Lorentz vacuum electrodynamic^ 

H^rF^M^F. (66) 



Remember that in Ricci calculus the excitation is defined according to TC- 
e'^"' Hki/2 . Then W^mc = - fV^F'\ see ^i. 
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7 Non-minimal coupling involving curvature, 
nonmetricity and torsion? 

7.1 Non-minimal coupling violating charge and/or flux conservation 

The Maxwell equation dH = J reflects (and comes from) the electric current 
conservation, dJ — 0, see Sec.|[ By modifying the left hand side of dH = J, 
one can arrive at a model violating charge conservation. Such a modification can 
typically originate from a non-minimal coupling of the electromagnetic to the 
gravitational field. Given the torsion 2-form T", one can consider, for example, 
the field equation 

dH + a (e„jT") AF + /3*(z?„ AT") AH = J, (67) 

or, with the nonmetricity 1-form Qap '.= —Dgap and the Wcyl 1-form Q := 

dH + jQ AH + S*{i}'^ A^'^ AQaf)) AH = J . (68) 

Similar non-minimal terms could emerge in dF — 0. However, curvature depen- 
dent terms cannot be accommodated at the level of the Maxwell equations, since 
the contraction of the indices produces always a form of even rank whereas the 
Maxwell equations are represented by 3- forms, i.e., by forms of odd rank. In any 
case, violating charge or flux conservation is not possible without giving up most 
of the experimentally established structure of the theory of electromagnetism. 
Therefore we will not follow this path. 

Incidentally, there are some papers in the literature in which the conventional 
vacuum constitutive law H ~ *F is uphold, but the Maxwell equations are 
coupled to torsion in an inconsistent way. A closer inspection of the papers 
[^Sjjl^,^ shows that the proposed "non-minimal" coupling of torsion to the 
electromagnetic field is void of physical contents. In fact, torsion drops out if the 
algebra is done correctly. 

Another procedure comes to mind if we talk about the violation of the conser- 
vation laws. Hojman et al. | ]36| introduced a new scalar field ip{x), the tlaplon, see 
also Mukku & Sayed [Q. The gradient dip of the tlaplon was put proportional 
to the trace part (2)t" := ,9" a T of the torsion T", see 0; here T e/jjT''. 
In fact, we have T — ^dip. 

Superficially, the axial scalar a(x) (Abelian axion) and the scalar tlaplon ip{x) 
may look similar. However, the axion already emerges from spacetime viewed 
as a differential manifold as soon as a linear constitutive law is assumed for 
electromagnetism, whereas the tlaplon can only be introduced if the differential 
manifold in equipped with a linear connection. In other words, the axion is a 
pre-metric and a pre-connection animal, the tlaplon, in contrast, needs to be 
'housed' in a linear connection. 

Moreover, as we saw, the axion respects the conservation laws (and pleases 
us thereby), whereas the tlaplon defies these rules and appears as an anti- 
electromagnetic creature. The electromagnetic field is not defined in the conven- 
tional way, namely by = dA. Instead, Hojman et al. define it via a "covariant" 
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derivative according to 

F^F~T"Aa^F+^T/\A = F+ l-dipAA. (69) 

Thus the "electromagnetic" Lagrangian becomes 

F A*F ^ F A*F + dip AAA*F AAA*{T AA) . (70) 

The last term does not contribute to the excitation H = —dL/dF, but the 
second term produces a contribution 

r-* {dip A A). (71) 



We can compare now the two contributions from the axion (48) and the 
tlaplon d?^). They are reminiscent of each other since one is equal to the Hodge 
dual of the other. Therefore, in the tlaplon case, besides a connection, we need 
additionally a metric. But the most decisive difference is, as can be read off from 
([toI), that the Maxwell equations get amended and axiom 1 and axiom 3 are no 
longer valid. 

7.2 "Admissible" non-minimal coupling 

The message is then that a change of the Maxwell equations dH — J , dF — 
is to be avoided, unless one allows for a violation of electric charge or magnetic 
flux conservation. By introducing the metric g^^ into the constitutive law, one 
gets a smooth and natural transition from special relativity to general relativity 
and to gauge theories of gravity. In the constitutive law for vacuum one could 
imagine, along with the contributions depending only on the metric, couplings 
like H 

+ A4 {Ric'^'^g^^^ ~ Ric^^'^g^^') + A5 R g'^^g^^^ + Aq R e*^'*^' (72) 

without violating the conservation laws. Here R := e" \ep\Ra^ is the curvature 
scalar, and we denote the right or Lie dual of an so(l, 3)-valued form by 

Choose, for example, -^^^''^ = 2y^—gfo{l + 0^ i?) g^I'^g'l^, then the inhomoge- 
neous Maxwell equation would read 

d*F + P^d {R*F) = ^J, (73) 
Jo 

i.e., a coupling of curvature and electromagnetic field strength would be possible. 
However, one had to introduce a new natural constant with the dimension of 
1/3] = 1/length. On the level of the Lagrangian, this coupling would be non- 
minimal, 

Kon-m = - y (1 + R) *F A F , (74) 
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but such an ansatz would not spoil the fundamental principles of electrodynam- 
ics; it would seem to be the most natural way of achieving an i? F-coupling. 
Goenner | |2^ , see also the literature given there, derived this non-minimal La- 
grangian from some fundamental principles, like the existence of a decent New- 
tonian limit. However, in his view, such a model violates charge conservation. 
We disagree with him on this point. 

We want to stress that one cannot achieve a similar non-minimal coupling to 
the torsion of spacetime. First of all, the Maxwell equations are independent 
of torsion. By means of the constitutive law one maps the 2-form F to the 2- 
form H. The curvature Ra^ is a 2-form of type (1, 1), i.e., it carries two Gi(4, R) 
indices, whereas the torsion 2-form carries only one index. Therefore, by con- 
traction, we cannot get a scalar out of the torsion. A coupling like (e^jT") A *F 
is not possible, since this is a 3-form. However, higher powers in T" would be 
possible such as 

ir = /o [l + 7'*(r"AT„)] A*F (75) 

or 

H = /o [1 + (54*((ejr") A iep\T^) A (e^jT^) A (e^jT*))] A *F . (76) 

Here [7] = [S] = 1/ length. Accordingly, it is not too difficult to introduce a 
coupling of torsion to electromagnetism. However, the price one has to pay is the 
introduction of new natural constants 7 and S. In other words, even if possible, 
we don't take such models too seriously. 

Also non-minimally coupled nonmetricity could be installed by additional 
quadratic pieces such as 

H^f, [l + ^'(e„Jg"^)(e/3JQ^)] *F . (77) 

Therefore, there are quite a number of different "admissible" options available 
as soon as we allow non-minimal couplings to arise. 

8 Outlook 

Using astronomical observations on the propagation of light, the upper bounds 
for non-minimal coupling effects should be determined in a systematic way, as 
is done, for example, by Haugan and Lammerzahl p^ . For such a purpose, we 
will develop the geometrical optics limit of the Maxwell equations dH ~ 
J, dF = and will use particular constitutive laws, as, e.g., the linear law. 
Possible couplings to curvature, torsion, and nonmetricity should come under 
sharper focus in this way. Non-linear effect a la Ayon-Beato & Garcia should be 
investigated in the context of, say, the metric-afRne gauge theory of gravity. 
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